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We study how universality classes of 0(TV)-symmetric models depend continuously on the di- 
mension d and the number of field components TV. We observe, from a renormalization group 
perspective, how the implications of the Mermin- Wagner-Hohenberg theorem set in as we gradually 
deform theory space towards d = 2. For fractal dimension in the range 2 < d < 3 we observe, for any 
TV > f , a finite family of multi-critical effective potentials of increasing order. Apart for the TV = 1 
case, these disappear in d = 2 consistently with the Mermin- Wagner-Hohenberg theorem. Finally, 
we study 0(TV = 0)-universality classes and find an infinite family of these in two dimensions. 
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Introduction. Our modern understanding of quantum 
or statistical field theory is based on the ideas put for- 
ward by K. Wilson and formalized within the framework 
of the renormalization group (RG) This approach 
considers all possible theories describing the quantum or 
statistical fluctuations of a given set of degrees of free- 
dom, the fields, subject only to the constraints imposed 
by symmetry and dimensionality. This defines what we 
call "theory space", the space of all possible theories. The 
process of quantization on one side, or averaging on the 
other, is then seen as a trajectory connecting the bare 
action or the Hamiltonian to the full quantum or statis- 
tical effective action. This trajectory can be of finite or 
infinite length (with respect to the RG time); in the first 
case one is performing an effective field theory calcula- 
tion, while in the second case one needs an ending point 
in theory space: this usually is a fixed-point. RG fixed- 
points describe scale invariant theories, i.e. theories that 
don't have a characteristic length, where fluctuations on 
all length scales are equally important. These theories are 
like lighthouses that permit us to shed light on the struc- 
ture of theory space. They attract or repel surrounding 
theories giving rise to universality, a phenomenon that 
underlies both non-perturbative renormalization and the 
understanding of continuous phase transitions. Once we 
find all fixed-points we can reconstruct the general (topo- 
logical) properties of the RG flow in theory space and 
acquire a deep understanding of that class of models. A 
paradigmatic example of this is the c-theorem Q, which 
describes the RG flow between two dimensional theories. 

Even if important information about two dimensional 
theories comes from exact results for particular lattice 
models, our ability to predict the universal features of 
two dimensional continuous phase transitions comes from 
the understanding of the RG structure of theory space. 
On the other hand, three dimensional systems are much 
more difficult to treat exactly; also here many analytical 
insights come from the RG study, while otherwise one has 
to resort to numerical methods. Deep insights, such as 
the role played by conformal symmetry in constraining 
statistical fluctuations, are also naturally embedded in 
the larger framework of RG analysis Q . 

In this Letter we show how another fundamental 
and broad result like the Mermin- Wagner-Hohenberg 
theorem [4j (see [5( for a review) , which states that there 
cannot be continuous phase transitions in d — 2 systems 
characterized by continuos symmetries, fits in the RG 



picture. We will do this by studying linear scalar O(N)- 
models, a class of theories that has many applications 
to low dimensional systems. Despite its relevance, there 
is no complete description of how universality classes of 
0(TV)-models depend continuously on both d and TV. In 
this Letter we give such a description by studying scaling 
solutions of the effective average action As a result 
we find many new TV > 2 universality classes describing 
multi-critical models in fractal dimension between two 
and three. In the TV = case we observe an infinite 
number of fixed-points in d = 2, analogue to the TV = I 
minimal models series. 

Flow equations. The effective average action (EAA) 
r^fyj] is a functional that depends on the infrared scale k 
and that interpolates smoothly between the bare action 
for k — > oo and the standard effective action for k — > @ . 
The EAA satisfies an exact RG equation [7] that describes 
its dependence upon changes of scale; this equation can 
be used to set up a framework where to concretely imple- 
ment the RG ideas discussed above. It is generally quite 
difficult to follow exactly the flow of the EAA through 
theory space and to find the relative fixed-point func- 
tionals: approximations are needed. One that retains all 
the important information about the structure of theory 
space is the one where all one-particle-irreducible (1PI) 
vertices of the EAA are evaluated at zero momenta. This 
defines the running effective potential Uk{p) which is a 
function of the 0(TV)-invariant p — \^p 2 ■ The running ef- 
fective potential can equivalently be defined by evaluating 
the EAA at a constant field configuration ip(x) = <p where 
Tk[<p] = ^Uk(p). In this approximation the RG theory 
space is projected onto the functional space of effective 
potentials. This space is still infinite dimensional and, at 
least at the qualitative level, all 0(TV)-universality classes 
of the full theory can be found by determining the relative 
scaling solutions (fixed-points). 

In terms of the dimensionless (running) effective po- 
tential U k {p) = k- d U k (p), with p = k-( d - 2+ ''tip, a scal- 
ing solution dtU*{p) = satisfies the following ordinary 
differential equation 
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Figure 1: Anomalous dimensions 772 as a function of d for 
(from above) JV = 1,2,3,4,5,10,100. In the inset we show 
the anomalous dimensions in the range 3 < d < 4. 

where c^ 1 = {Air) d / 2 T(d/2 + 1). The anomalous dimen- 
sion r\ fixes the scaling properties of the field at a particu- 
lar fixed-point; to lowest order the value of the anomalous 
dimension is related to the dimensionless effective poten- 
tial by the following relation 

V = c d - 72 » I 2 ) 

1 + 2 / 5 l>( / 5 ) 

where po is the (absolute) minimum of the dimensionless 
effective potential, i.e. f7»(/5o) = 0. 

Every scaling solution, together with its domain of 
attraction, represents a different universality class; thus 
by finding the solutions of the system composed of ([1} 
and © one can determine all 0(JV)-universality classes. 
Differently from other implementations of RG ideas, all 
of the analysis can be made leaving d and JV as free 
parameters, permitting us to study how theory space 
depends on these. 

Mermin- Wagner- Hohenberg theorem. We solve the 
fixed-point equations ([J} and © by the iterative method 
proposed in |3j- Given d,N and an initial ansatz for the 
anomalous dimension one can numerically integrate equa- 
tion ([T]), subject to the initial conditions Ui(0) = and 

U* (0) = d(d+2)^i+cr) ' ^° ^ nc ^ ^hat f° r mos t values of the 
parameter a the solution will end up in a singularity at a 
finite value of the dimensionless field; let's call this value 
Ps' N ' v (<r). Requiring a scaling solution to be well defined 
for any p £ K restricts the admissible initial values of a 
to a discrete set {o-^'f' v }. As we will see, all these solu- 
tions correspond to multi-critical potentials of increasing 
order; we will use the subscript i to specify the number 
of minima of a scaling potential. One can construct a nu- 
merical plot of the function pf N,v (a) to find the a^' n as 
those values where it diverges or has a "spike". Then one 
inserts the solution of the ODE ((1) obtained from a value 
o-^^ ,v , and thus corresponding to a i-critical potential, 
in (|2J| to find the next value for the anomalous dimension, 
and iterates the procedure until convergence. In this way 
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Figure 2: 0[d, N) = rj 2 (d, N)/r) 2 (d,l) as a function of N for 
(from above) d — 2.1,2.05,2.01,2. The anomalous dimen- 
sion can be interpreted as the order parameter of a Mermin- 
Wagner-Hohenberg continuos "phase transition". The con- 
trol parameter of the transition is JV and the critical value 
is JV, = 2. 



it is possible to construct the function r\i (d, JV) describing 
how the i-th scaling solution depends on the dimension 
d and on the number of components JV. 

By this technique we can follow the evolution of the 
RG fixed-points as we vary d and JV. For d > 4 and for 
any JV we find only the Gaussian fixed-point; at d = 4, 
the upper critical dimension for the 0(JV)-models, the 
Wilson-Fisher fixed-point starts to branch away from the 
Gaussian fixed-point. In d — 3 these fixed-points de- 
scribe the known universality classes of the Ising, XY, 
Heisenberg and other models. Approaching d = 2 one 
clearly observes that only the N = 1 anomalous dimen- 
sion continues to grow: for all other values of N > 2 the 
anomalous dimension bends downward to become zero 
exactly when d = 2. This is a non-trivial fact, not ev- 
ident from the structure of equation ([I]), telling us that 
only the 0(AT)-model with discrete symmetry (JV = 1) 
can have a second-order phase transition in two dimen- 
sions, while all the 0(iV)-models with continuos sym- 
metry (JV > 2) cannot. This result, that here emerges 
from the RG analysis alone, is commonly known as the 
Mermin- Wagner-Hohenberg (MWH) theorem [4|. In Fig- 
ure 1 we plotted the anomalous dimension as a function 
of d for the particular values JV = 1, 2, 3, 4, 5, 10, 100; this 
figure clearly represents the way in which the MWH man- 
ifests in the RG framework. This analysis can be seen as 
a RG confirmation of this important result and can be 
the starting point for a new rigorous proof of it. Note 
also that, as expected from the exact solution [9|, the 
anomalous dimensions tend to zero for JV — > 00. 

That there is no continuos phase transition in d = 2 
for the JV > 2 models can be confirmed by the analysis of 
the temperature critical exponent ^(d, JV). This allows 
us to distinguish the Spherical model, related to the JV — > 
00 limit, from the Gaussian model, both having 77 = 0. 
We find Q that only the JV = 1 model has a finite 
in two dimensions, in all other models it diverges upon 
approaching d = 2, as is found in the JV — > 00 case where 




one knows exactly that ^(d, oo) = — yrg- An RG study 
by EAA techniques of this last case can be found in Q . 

The critical case TV = 2 is known to have a distin- 
guished behavior (loj . In this case one can observe all 
the distinctive properties of the Kosterliz-Thouless phase 
transition by studying the properties of the RG flow; we 
refer to [ll| for more details. 

It is worth noting that our result for the anomalous 
dimension 772(2,1) = 0.234 is in good agreement with 
the exact result r\ = 0.25, but we must mention that to 
provide an error estimate on this value one needs to go 
to higher orders of the derivative expansion [l2| ■ 

To better discriminate between theories which can un- 
dergo a continuous phase transition in c?* = 2 and those 
which cannot, we extend the analysis of scaling solu- 
tions to non- integer TV; in particular we want to see what 
happens around the critical value TV* — 2. The MWH 
theorem tells us that at d = d* the "order parameter" 
0(d,N) = i] 2 (d,N)/r]2(d, 1) is zero for TV > TV, and non- 
zero for TV < TV*; but it tells us nothing about the "order" 
of the transition. For example, is the transition "contin- 
uous" in TV or not? Figure 2 shows that the RG analysis 
can say a lot more about this. First, we see that 0(d, TV) 
evolves continuously with TV across TV*; second, we see 
that 0{d, TV) can be written in a scaling form around the 
transition point (d*,TV*) = (2,2); in particular we can 
write the following scaling relation: 



N.-N 
N* 





TV 
TV 



N- 
N+ 



(3) 



where we introduced a new scaling exponent 0. A fit from 
the data displayed in Figure 2 gives the estimate G « 0.98 
which is quite close to one. Relation ^ tells us how 
theory space deforms as we vary the external parameter 
TV. An interesting question might be if relation (|3j) is 
universal, in the sense that the value of O is independent 
of the details of the implementation of the RG procedure 
but rather describes an inner property of the set of theory 
spaces parametrized by TV. One can make a reasoning 
similar to the above by keeping TV fixed at its critical 
value but varying d around <i*: 
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In (HJ) we included the information, taken from [T^|, that 
77 remains zero for TV > TV* and d < d*. A fit from the 
data displayed in Figure 1 gives the approximate value 
A ss 1.86. Finally, we remark that our analysis works 
when the first term on the rhs of ([1} is non-zero, i.e. 
when r\ > 2 — d. This fact prevents us from performing a 
complete analysis in the range 1 < d < 2, where indeed 
studies of 0(TV)-models on fractals have shown that the 
MWH is still valid [lj|. 

Multi-critical 0(N)-models in fractal dimension. 
When new universality classes appear by branching from 
the Gaussian fixed-point it is easy to determine the 
relative critical dimensions, since the argument based 
on canonical dimensions is valid. In particular, the 
i-th multi-critical scaling solution appears at the crit- 
ical dimension d Ct i = 2 + jzj] for integer values of 
i = 1,2,3,4,... we find d c<i = 00, 4, 3, §, f, f , |, f , |, ... 
Each of these represents the upper critical dimension of 
the i-th multi-critical universality class. At these dimen- 
sions we see non-trivial fixed-points branching from the 
Gaussian even when TV > 2. These fixed-points indeed 
correspond to multi-critical potentials of increasing or- 
der, i.e. with i minima when expressed in terms of the 
variable 

The critical dimensions d c ^ accumulate at d = 2 
and thus one may naively expect to find, for any TV, 
infinitely many universality classes in this dimension. 
Our analysis shows instead, see Figure 3 for the cases 
i = 3,4,5 and TV = 1,2,3,4, that this happens only in 
the TV = 1 case, where the multi-critical fixed-points 
approach, in the limit d — > 2, the fixed-points repre- 
senting the minimal- models found in CFT Q. For any 
other TV > 2 we find that, consistently with the MWH 
theorem, the multi-critical scaling solutions, present 
in the range 2 < d < 3, are instead absent in d = 2. 
This fact is a strong check of the general validity of 
the MWH theorem, which our analysis indicates is also 
applicable to multi-critical phase transitions. On the 
other side, we predict the existence of a whole family of 
0(TV)-universality classes in fractal dimensions between 
two and three. To our knowledge these universality 
classes are new. 

The TV — > limit. In this last paragraph we study the 
TV — s- limit, that describes, through the correspondence 




Figure 4: r/2 as a function of N is continuos both in d = 2 
(upper curve) and in d = 3 (lower curve). The limit N — > is 
well defined and both curves approach zero at N = —2, where 
the theory is exactly solvable with n = 0. 
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Figure 5: Anomalous dimensions r\i as a function of d for the 
first four N — multi-critical scaling solutions, i.e. for (from 
above) i = 2, 3, 4, 5. 



found by De Gennes fl4| . the universality class of the 
self-avoiding random walk (SAW). Figure 4 shows the 
anomalous dimension as function of N in the interval 
between —2 < N < 2.5 for the cases d = 2 and d — 
3. The figure shows that the anomalous dimensions are 
continuous in the whole range; this is an indication that 
the N — > limit is continuous, and thus well defined. 
Figure 4 also shows, interestingly, that both the d = 2 
and d = 3 curves tend to zero as N — > — 2 where indeed 
the model is exactly solvable in both dimensions and has 
vanishing anomalous dimension 

As for the other integer values of N, we find multi- 
critical scaling solutions for N = 0. The interesting 
thing here is that these solutions survive in infinite 
number when d approaches two. A plot of the first 
four anomalous dimensions is shown in Figure 5; these 
are numerically very similar to those of the N = 1 
models (which can be found in Figure 1 and 3). This 
similarity is expected, as one may see by inspection of 



Figure 4. Even if the anomalous dimension is not a 
relevant physical parameter in the SAW correspondence, 
we can use scaling relations to relate it to the physical 
critical exponents v and 7 in order to make some 
checks. In d = 2 one finds the exact values v = | and 
and so r\ = 2 - J = ^ ~ 0.208; we find 
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772(2,0) = 0.232. In d = 3 one finds from Monte Carlo 
simulations the values v = 0.587 and 7 = 1.157 [l||, 
and so T] = 2 - 2 ~ 0.029; we find r/ 2 (3, 0) = 0.04. As 
we said before, we cannot extend our method to d < 2 
to compare with exact SAW critical exponents found 
on fractals [17]. In any case, our analysis suggests that 
there is a countable family of 0(N = 0) -universality 
classes in two dimensions. To our knowledge these are 
novel and may describe multi-critical phase transitions 
of some polymeric system. One may wonder if they are 
related to any CFT. 

Discussion and outlook. In this Letter we studied how 
universality classes of linear scalar theories with O(N)- 
symmetry vary continuously with the dimension d and 
with the number of field components N. We followed the 
evolution of RG fixed-points, represented here by scaling 
solutions of the RG equation for the effective potential 
((1}, through theory space as we varied these parameters. 
As in |H, even if all our analysis was based on the study 
of a simple ODE, we were able to observe a very rich 
behavior. 

Above four dimensions and for any N we found, as 
expected, only the Gaussian universality class, while at 
d = 4 we observed the Ising (or Wilson- Fisher) universal- 
ity classes appear. In fractal dimension between two and 
three we found non-trivial fixed-points for all N: these 
are novel universality classes that in principle can be ob- 
served both in theoretical models on fractal lattices or in 
real physical systems. 

Approaching two dimensions we observed the RG man- 
ifestation of the Mermin- Wagner- Hohenberg (MWH) the- 
orem: only the N = 1 universality classes survived down 
to d = 2, while all the N > 2 ones disappeared. By con- 
sidering (d, N) as real parameters near (2, 2) we found 
that the transition described by the MWH theorem, be- 
tween theories that can undergo a continuous phase tran- 
sition and theories that cannot, is continuous and that the 
anomalous dimension, which can be seen as analogous to 
the order parameter, can be written in scaling form at 
the critical point (2, 2). Our analysis revealed how differ- 
ent theory spaces parametrized by N are related to each 
other; this information gives a deep RG understanding of 
the MWH theorem and can be used as the starting point 
for an extension of it. 

Finally, we studied the N — > limit; we found that 
it is continuous around = and we gave evidence of 
new 0(N = 0) -universality classes in d = 2. These are 
analogous to the CFT minimal models series and may 
describe particular multi-critical transitions of polymeric 
systems. 



[1] K.G. Wilson, Rev. Mod. Phys. 47 (1975) 773; K.G. Wil- 
son, Rev. Mod. Phys. 55 (1983) 583. 
[2] A.B. Zamolodchikov, JETP Lett. 43 (1986) 730, [Pisma 

Zh. Eksp. Teor. Fiz. 43 (198 6) 565]. 
[3] A. Codello. larXiv: 1204.38771 [hep-th], 
[4] N. D. Mermin and H. Wagner . Phys. Rev. Lett. 17(22) 

(1966) 1133; P. C. Hohenberg, Phys. Rev. 158(2) (1967); 

Coleman, S. (1973) Comm. Math. Phys. 264 (30819) 259. 
[5] A. Ge lfert and W. Nolting , J. Phys. Con. Mat. 13 (2001) 

R505, |cond-mat/0106090| 
[6] J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363 

(2002) 223, |hep-ph/0005122| 
[7] C. Wetterich, Phys. Lett. B 301 (1993) 90. 
[8] A. Codello and G. D'Odorico, in preparation. 
[9] N. Tetradis and D.F.Litim, Nucl. Phys. B 464 (1996) 492, 

|hep-th/95 12073] ; T.R. Morris and M.D. Turner, Nucl. 

Phys. B 509 (1998) 637, [hep-th/9704202] . 



[10] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6 (1973) 
1181. 

[11] M. G rater and C. Wet terich, Phys. Rev. Lett. 75 (1995) 
378, [hep-ph /9409459|; G. Von Gersd orff and C. Wet - 
terich, Phys. Rev. B 64 (2001) 054513, [hep-th/0008114] . 

[12] L. Canet, B. Delamotte, D. Mouhanna and J. Vidal, 
Phys. Rev. D 67 (2003) 065004, |hep-th/0211055] L. 
Canet, B. Delamotte, D. Mouhanna and J. Vidal, Phys. 
Rev. B 68 (2003) 064421, [hep-th/0302227l 

[13] D. Cassi, Phys. Rev. Lett. 68 (1992) 3631. 

[14] De Gennes, Phys. Lett. A 38 (1972) 339. 

[15] B. Nienhuis, J. Stat. Phys. 34 (1984) 731. 

[16] A. Pelissetto and E. Vicari, Phys. Rep. 368 (2002) 549. 

[17] S. Elezovic, M. Knezevic and S. Milosevic, J. Phys. A 20 
(1987) 1215. 



